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THEOHETICAL SOLUTION OP PROFILE DRAG* 
By J. Pretsch 



SUMMARY 



After a survey of the customary procedures for ap- 
praising the profile drag in which pressure drag was dis- 
counted and of the methods for computing the laminar and 
turhulent friction flow, the author proposes a method hy 
which the pressure drag can he computed with the a.id of 
the displacement thickness of the frictional layer. The 
method is restricted to the case where the effects, caused 
hy separation of frictional layer, are small. Then the 
t tal profile drag can he expressed solely hy quantities 
derived from the velocity distrihution in the frictional 
la.yer immediately at the trailing edge. It is merel>*' 
assumed therehy that the mixing losses originating over 
the short length-in the wake are negligihle until the 
pressure reaches its end value. The proposed method is 
applied to seven symmetrical Karman— Tr ef f t z profiles at 
zero lift and varying position of transitional region. 
The actual positi.on of the trans it i ona! r egi on is deduced 
hy comparison with the measured drag coefficients. Judged 
hy these mathematical results and the availahle test data 
the inference is drawn that the position of the transitional 
region is principally dependent upon the Reynolds numher 
Re,5 referred to momentum thickness, hut for the rest al- 
most independent of the pressure gradient while being ; 
materially affected hy the degree of turbulence and sur- 
face roughness; the pressure gradient becomes naturally 
effective in the quantity Re^ itself. This supposition 
is utilized to predict the transit.i'onal' "region and hence 
the profile drag coefficient of a smooth wing at moderate 
degree of turbulence and very large Reynolds numbers 
(lO'^ < Re < 10 ) and to compute the friction and pressure 
drag. It seems as if the percent proportion of. the pres- 
sure drag to the total profile drag increases with the 
Reynolds number. 



*"Zur the ore t is Chen Berechnung des Pr of ilwiderstandes , " 
Reprint from Jahrbuch 1938 der deutsche Luf tf ahrtf or s ch- 
ung , pp. I 60 — I 81. 
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In c one lus i'on. the most urgent .prqt lems of the future 
are enumerated, with the solution of which an improvensent 
in the theoretical profile drag computation may be antic- 
ipated. . . . , . . 

Notation 

Velocities (m/ s ) : 

uoo flying speed 

velocity at the outer "boundary of the frictional 
layer 

u tangential component of the velocity within the 
frictional layer 

V normal component of the velocity within the fric- 
tional layer 

Lengths (m) 

X coordinate of a point of the wing in flight direc- 
tion measured from the stagnation point 

y coordinate of a point of the wing in span direction 

z coordinate of a point of the wing at right angles to 
X, y plane 

s arc length of a point of the wing measured from 
stagnation point 

n distance from surface of wing 

t wing chord ■ 

■^max "Maximum profile thickness 

"b span 

6. Pohlhausen's boundary layer thickness 

5* displacement thickness 

* momentum thickness 
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Angles {. .- 

P angle between profile tangent and flight direction 
cp trailing edge . . 

Pressure (kg/m^) 

p static pressure 

q dynamic pressure 

t = cLt\ difference in dynamic pressures at the friction 
layer "boundary and at the point nf momentum thick- 
ness 

Tq shearing stress 

J'orces per unit length of span (kg/m) 
W force in flight direction 
Wp profile drag 
Wj. frictional drag 

pressure drag 

Other physical quantities: 

p air density j i 

t_ m* J 



M. air viscosity ^£.1 
u kinematic viscosity i | 



E yield of 
2 - IT Cs) 



a source ■^ 



m 



3-1 



s 
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Bimens ionless quantities: 

d maximum thickness in ratio . to .wing chord 
f maximum camber in ratio to ving chord .-- 
c^ lift coefficient 



coefficient of profile drag 



c^^ coefficient of friction drag 
c,,,^ coefficient of pressure drag 

■Uoo t 



Re = 



He 



6* 



He 



. V 

V 

uo <^ 



Reynolds numbers 



1? 



* u.* 6* 
Re„* = -2 



X = ^ d_uo 
~ V d s 

n = 1 



* _ 



L 



Pohlhausen's form parameter 
s 

G-ruschwitz' form parameter 



Uq d s 



L ^ 



^0 S 



P Uo' 



V 



/4 

Uikuradse 's form parameter 
dimensionless shearing stress 



Indices : 

u transitional region of friction layer 
a separation point of friction layer 
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I. IWTEODUCTIOIT 
Division of Profile Drag in Friction and Pressure Drag 

"She aim of the present inye s t igat i on is to compute 
the profile drag of a wing or of any other cylindrical 
"body with streamlined profile without resorting to meas- 
urements. There is no or only ins ignif'icant separatlo:n 
of the friction layer on such profile forms at small an- 
gles of attack. The profile drag is pr edomina;ntly s.kin 
friction and relatively small in c6mparison to "bodies 
with c ons iderah le . spread of separation^ The study is 
confined to such "hodies of low drag," and of these to 
such with approximately two— dimens ional flow. 

Since on the modern high— speed aircraft the profile 
drag can amount to atjout half of the total drag the 

problem of theoretical solution is accorded great prac- 
tical significance. It is rendered difficult hy the 
multiplicity of effects which are ahle to modify the mag- 
nitude of the profile drag. 

Whereas the induced drag depends solely on the de- 
sign of the wing at full scale, that is, its contour and 
its profile and is computable from the drawing heard 
design itself for all flight cases, the profile drag is 
also dependent upon the method of construction, such as 
r oughnes s, wavines-s in covering, position of rivet head 
rows, flap gaps, etc., and in flight close -to the ground 
on local weather conditions, that is, the degree of tur- 
bulence 'of the air streaming past the wing. 

In order to compute the profile drag the- wing must, 
in the true sense of the word, be put under the magnify- 
ing glass, by reason of the fact that the drag is largely 
induced by the processes taking place in the iII^mediate 
pro3cimity of the wing surface, namely, in the- friction 
layer . 

The air, by adhering to the surface of the wing be-, 
cause of its viscosity, transmits a shearing stress T.^ 
to it; this stress acts, in two-dimensional flow-, in 
direction of the tangent to the wing profile. 

If, starting at the stagnation point, the arc length 
is denoted with s and the coordinate of a point of the 
wing in flow direction with x, the resultant Wj. of 
all shearing stresses in flight direction referred to unit 
length of span b (fig. 1) is: 
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Vj. =: jT T.^ cos 3 ds = . dx (1) 

1+ II 1+ II 

where ^ denotes the, angle between profile tangent and, 
flight direct.ioni-.and suT3scripts I and II indicate 
that, the integration from the stagnation point is to be 
extended once over the suction side and once over the 
pressure side. The force Wp h is termed' the friction 
drag of the wing; if represents the first portion of the 
profile drag. 

The second loss inducing effect' of the air viscosity 
is the changed pressure dis trihut i on, along the total wing 
contour due to the pushing aside of the potential flow. 
These losses occur even on low drag bodies, hut hecoffie 
especially apparent when the friction layer separates, 
because it no longer possesses sufficient kinetic energy 
to meet a strong pressure rise of the potential flow. 
Behind the separation point th6 pressure actually pre- 
vailing does not rise again appreciably; therefore a low- 
pressure region builds up, and the wing is sucked backward 
,ag-ainst the direction of flight. 

With p denoting the static pressure and z the 
coordinate of . a profile point at right angle to the direc- 
tion of flight, the- resultant 6f all normally acting 
forces in flight . direct ion , referred t« unit length of 
span b (fig. 2) is : 

= r p sin 3 ds = r p dz (2). 
II 1+ M 

the subscripts I and II having the same meaning. as in 
equation (1). The force ¥^ b is termedthe pressure 
drag .of the wing. 

The shearing stress "^o' through which the frictional 
drag is computed according t& (l) can at least be approxi- 
mated by the calculation of the friction layer. And, if it 
succeeds in reducing the pressure drag to quantities that 
characterize the friction layer, the problem of .profile 
drag reduces to that of defining the friction layer along 
the wing. This method is to be developed in the present 
report. 
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II. -PEEVIOUS THEOREMS POH PEOyiLE D£A& APPEAISAL 
Seduction to Profile Drag of Plat Plate 



The only "body on which the profile drag has "been 
amenahle'to solution is the thija , flat plate tangent ially 
exposed to a flow of constant spSeed Uos, . Its^ pressure 
drag is zero because all pressures are at right angles 
to the direction of motion; its profile drag consists 
therefore exclusively of frictional drag. To compute 
this frictional drag, three cases must he differentiated 
theoretically, depending upon the character of flow within 
the friction layer. 

a) The flow- within the friction layer past the entire 
plate is laminar. In this in-stance the drag coefficient 
of the plate referred to the "base surface is, according 
to Blasius (reference l): 

= ;. Re = - ' (3) 

where t is plate length; i- , kinematic riscosity. This 
law holds true up to ahout Re = 5 x 10^ . 

h ) The friction layer along the plate is turbulent. 
In this case 



'W, 



0.91 (4) 



(logio ^e) 



3.58 



up to any large values of Re, according to Schlichting 
(reference 1 ) . 

c) The flow is first laminar, then turbulent. The 
transition from one to the other takes place over a dis- 
tance of finite length, hut may he visualized for the drag 
calculation as being localized in a point whose distance . 
X from the leading edge can be defined by 

^ X 5 ■ 

~ 5 X 10 (5) 



V 



whereby 



0.91 .3400 - 



(logio Re)^'^^ 
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likewise .bjolds for any Reynolds number, according' fo 
Prajidtl (referen.ce 2). „. 

The simplifying assumption of sudden reversal signi- 
fies that on the transitional distance the actual shearing 
stress- distrihuti on ("fig. 3, solid curve) is replaced, by 
the laws.foT. purely laminar or purely turbulent friction 
layer flow (fig. 3., -dot'ted curve), which has no effect on 
the integral . value of ; the drag. 

The formulas cited refer to the smooth plate; although 
those for rough plates have aliso been studied exhaustively 
by Prandtl and Schlicht ing ' (reference 3 ) . Since, the pro- 
file drag law of the thin, flat plate of ' the "friction 
sheet" as we shall call it for short, is known, it is ex- 
plicable why it has been ' repeat edly attempted to estimate, 
the prof ile ■ drag- of a wing, even that of an airship hull > 
by substitution of ah "equivalent friction sheet" (Jones, 
Hoerner , Betz, Bock,Dryden, and Kuethe); 

The .roughest approximation is obtained iirith Jones 
(r ef ere nc e • 4 ) when the wing is replaced by a just— as— fast 
flying smobth plate of the same depth, the base of which 
is equal tb half the total wing surface. 

Hoerner (reference 5) computes the profile' drag of 
modern, high-speed aircraft as friction drag of likewise 
identically fast but also just-as-rough friction sheets, 
conformal to the Prandt 1-S chl icht ing formulas for rough, 
flat plates (reference 3). 

Hoerner and Jones' approximations make no allowance 
for the relationship between shearing stress and pressure 
distribution, which Betz (reference 6) had attempted to 
estimate back in 1915, For, on assuming that the shear- 
ing stress is proportional to the square of the speed 
Ug on the outer boundary of the friction layer, the fric- 
tion drag- coefficient of the wing in ratio to that of the • 
equivalent sheet is: 

'^^r wing ■ ; / u„\^ <? 

i = / ( ^0 ) d - (7) 

Cw plate J \ ■n<=°/ t 

where t is wing chord; u^, the flying speed. 

An estimation of wing friction drag similar to that 
by Betz is afforded in the proposal by Beck (reference 7), 
namely, to take the f lying spee^d of .the equivalent fric— 
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tion sheet ■ eci.ua 1 to the "mean f Ipw vel oc i ty" .pas t -the 
wing, hut fails to g'ive general directions regarding the 
type of averaging. - 

Dryden and Kuethe .(ref erejice . 8) go.a..step farther 
in the appraisal of the profile drag of airship hulls, 
vrhich is mentioned here "because . t heir , me thod would, he 
equally applicahle to airfoils.: T.he area ...of the "substi- 
tute plate" -is computed first, its .width .,at ..'a give'n dis.— . 
tance from the front edge heing equal to^half the . oir.cum— 
ference of the airship ring at the same distance from the 
nose. This "flattened airship" is cut a.long that "ring" 
at which, the reversal from.laminar to. tllrbulent .. f rict ion 
layer flow is assumed. Oh the forehody of the plate the 
friction layer is laminar, on t he . af t erh .o dy^, turbulent. 
The profile drag of the "stern plate" Dryden and Kuethe. 
compute on the assumption that the f low. veloc ity every- 
where is equal to .the speed u,jo of the airship. But the 
profile-drag of the "nose plate" they compute on the as- 
sumption that the speed outside of the friction layer 
changes from one. ring to the next exactly as the velocity. 
Uq on the airship. Then a linear theorem for .the veloc- 
ity distribution within the laminar friction layer enables 
them to compute the shearing stress Tq and hence the 
drag of the forward portion of the friction sheet on. the. 
basis of Von Karman's momentum equation. 

This method of friction sheet analogy is now subject- 
ed to a critical analysis: 

a) One obvious defect is that the pressure drag can- 

not be solved at all .by theoretical computa- 
tion. This defect is so much more palpable as', 
the profile is thicker. 

But even the friction drag itself cannot be properly 
determined by this method and for the following reasons-: 

b) The shearing stress in the laminar as in the tur- 

bulent zone of the friction layer is definitely 
dependent upon the pr e s sure . di s t r ibut i cn along 
the profile contour. Dryden. and Kuethe allowed 
for this i.n the laminar port ion-,- but net. in the 
turbulent . ' ' '. 

c) The position of the t rans i t i onal. r.egi on depends, 

according to (5), not only on the Reynolds num- 
ber Re, but also on the pressure variation 
and on the degree of turbulence. 
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d) 



The angle p "between arc element vand flight di- 
rection is other than zferoon the airfoil; 
hence the shearing stress is effective with 
only one component, though on slender profiles 
■ the error involved is quite small-. 



However, the discrepancies enumerated of the charac- 
teristics of an airfoil from those of an equivalent flat 
plate may "be so pronounced as to make the application of. 
the analogy "between the two very restricted. 



Statistical Appraisal of Experimental Results 

This method consists of evaluat ing . exper imentally 
secured, values on a mult itude of airfoils statistically 
so as' to afford a practical empirical formula for compu 
ihg the profile drag for given profile parameters, if 
possi"ble in respect to angle of attack. In addition to 
this, a couple of tables from which eventual increments 
for degree of tur"bulence , surface condition, etc., could 
"be read off, would "be- des ira'ble . 

One pro"blem found to "be very deterrent in this method 
concerned the extent to whdch model tests are at all com— 
para"ble when carried out in different wind :tunnels , that 
is, different tur"bulence factors, different surface finish, 
and possi"bly different wing tip design (references 9, 10). 

Moreover, this method will not afford any insight 
into the mechanism of profile drag, perhaps in connection 
with friction and pressure drag, in role of the transi- 
tional region or the area of separation. Statistical cor- 
relation of recorded total drag values will simply lead to 
statements concerning the total profile drag, unless addi- 
tional measur ement s , as, for instance, of the pressure dis— 
tri'bution along the profile contour or of .the velocity 
distr i"but ion in the friction layer are available for each 
airfoil . 

Munk (reference 11) evidently in continuation of 
earlier studies (reference 12) gives the following rough 
rule for the profile drag coefficient 




(8) 



where d and f denote thickness a^d camber in ratio to 



NACA Technical Memorandum No. ,1009 



11 



wing chord, Ca the lift^ coefficient . • But never 

to he less than 0.006; this minimum value ' def ines the drag 
of a plate at Re = 10 , . 

Another formula "by Jacohs, Ward, and Pinkerton (ref- 
erence 13.) ohtained from .measurements . on 78 .different pro- 
files reads : 

' % .= °wp „,in ^ >wp' • . ; -(9) 

for symmetrical profiles (suhscript s): 



c 



v, , = 0.0056 + 0.01 d + 0.1 d^ (lO) 

"p s min 



where d is thickness in ratio to wing chord. 

For camhered profiles of the same thickness: 

c„ .. = ■ c„ ; + k;. k > 0 (11 ) 

*p min Wp s mm 

k is shown in the. diagrams plotted against position and 
amount of maximum camber. The angle of attack increment 
A c„ can also he read off from charts. 
P 

According to Doetsch and Kramer (reference 10) the 
measurements in the KACA var iah le— dens i ty tunnel, evaluated 
in reference (13), yielded, "because of the high turbulence 
factor (critical sphere characteristic Refc = 1-2 x lO^O 
and the "blunt wing tips involving an additional pressure 
drag at the tips, profile drag coefficients hy about 18 
to 35 percent too high. So the formula (9) should afiford 
at best a q.ualitat ively correct picture. 

Very comprehensive is the comparison of profile dra-g 
coefficients carried out by Glass (reference 14)-. He 
studied the change of minimum profile drag with the 
Reynolds number but disclaimed the angle of attack rela- 
tion. The effect due to turbulence he ttismisses .as' 
negligible on the basis of a report by Hoerner (Teferenice 
15). This is obviously due to some misinterpretation. 
Glass also ultimately reduced the profile drag of a wing 
to the drag of an equivalent, identically rough friction 
sheet, just as Hoerner did, but, contrary to Hoerner , he 
allowed for the effect of profile parameter by a correc- 
tion factor. 
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Glass first divided .the coefficient of the profile 
drag into two parts: 

% min = plate * ^ % ^l^) 

The first part represents the drag of the ecLui'valent 
friction sheet with the plan form and surface condition 
of the** wing. The second part is a function of the geo- 
metric determinative pieces of the profile which is again 
split into 

^ ?5wp = ^ °Wj. + Cwd • (13.) 

Here A Cvp denotes the cumulative friction drag "due 

to surface curvature," that is, the exce.ss of wing fric- 
tion drag relative to friction drag of equivalent plate; 
c^,^ is the coefficient of pressure drag. . 

From- the evaluation of extens ive .test data from 
G-erman, English, U.S., and Russian wind tunnels, G-lass 
established the following remarkah le- f ac t s : 



X8i 



a) Independent of Reynolds numher and airfoil ,shape 



const = 0.85 (14) 



This law is deduced from tests "by Fage , S'alkner and 
Walker (reference 16) on seven symmetrical Zarman— Tref f t z 
profiles, in which the pressure distribution was also 
recorded experimentally. A detailed discussion on this 
subject follows. • < 

"b ) Thickness d and camber f , both referred to 
wing chord, can be combined in a s ingle • ge ome tr ic param- 
eter, the s 0— called equivalent thickness" d* 

d* = d + 0.17 f2 (15) 

then it is found that by given. Re, the profile drag por- 
tion A c^^ grows r-inearly with the equivalent thick- 
ness d* ^ ■ . . • 

.. Ac^ ; = Z "(Re) -d*.. ■ " ■ (16) 

■ . P . ■ • . 



Division of the factor K in the product 



'NACA Technical Memorandum No. .1009 



13 



K (Ee) = (Re) ^ ' (Ee) (17) 

. Opiate . ■ 



gives , f inally ; 



P min 



= °w plate [1 + - (18) 



where c„ , . can "be computed "by the Prandtl— Schlichting 
Wplate 

formulas, n is represented as explicit function of Ee . 

G-lass 'quotes 3 x 10^ < Re < 9 x 10^ as the range of 
validity of his - empir ical formula (18); it is the very 
zone in which, on the smo o t h f r i c t ion sheet, the flow 
changes from laminar to turhulent . So the initial curve 
d* = 0 is exactly Prandtl's transition curve (6): 

But even G-lass ' formula, as meritorious as it is, 
can, "because of the discounted effects of turhulencS and 
limitation to minimum drag value, neither make claim to 
the desired perfection nor, because of the defective com- 
parability of the test values, even to satisfactory accu- 
racy. 

I.n the following it is attempted to compute with 
the means available the profile drag in theoretical man- 
ner, by' reducing its two proportions, the pressure and 
friction drag to quantities that characterize the fric- 
tion layer. Even though the computation of the friction 
layer which becomes necessary herewith, and that of the 
profile drag is in many points still beset with uncer- 
tainties, it nevertheless affords a valuable explanation 
of the question what circumstances are primarily deci- 
sive for the drag and so point out the possibilities of 
lowering the drag. 



III. FRICTION LAYERS ON A WINS 



First we secure a qualitative picture of the flow 
in the friction layer of a wing. Near the stagnation 
point the air stream will, be laminar. Nejct follows a 
transitional region in which the laminar flow turns tur- 
bulent, then comes a distance with fully developed tur— • 
bulence within, the friction layer, and, lastly, the po- 
tential separation, if the pressure rise is large enough, 
followed by a dead— air region extending from separation 
point to trailing edge. 
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Restrict ed' to "b odieg • of -• l-ow drag the problem of com- 
puting the friction layer resolves itself into ff^ur 
stages: ■ the laminar layer, the reversal point, the tur- 
bulent layer, and the p«int of separation. The solution 
is predicated first on pressure distrihution across the 
body surface. 



Pressure Distribution 

In order., to solve the profile drag of a wing without 
resorting t o" measurements it must be possible also to com- 
pute the pr es sure ■ dis tribut i on as. it really •ccurs under 
the effect of "air friction. This, .is • not. always possible 
to be obtained with desired accuracy. To be sure, there 
are methods! to determine the potential theoretical pressure 
distribution ,on .families of , pr of iles ( Joukowski » Karman— • 
Trefftz, Betz— Keune (reference 16)) and even on arbitrarily 
shaped profiles, (iref er ences 17 , ■ 18.) , ■ but the measured pres- 
sure distribution are not in . sat isf act ory agreement with 
those obtained by potential theory, especially at higher 
angles of attack. On the other' hand, Betz has pointed out 
as far back as 1915 (reference 6) that a much better agree- 
ment could be achieved if the circulation is secured from 
the measured lift rather than from the conditir^n of smooth 
trailing edge flov off. Admittedly, it results in flow 
around the trailing edge and hence in excessive low pres- 
sures. This obstacle has been removed recently by Pinker- 
ton (reference 19) while retaining Betz ' artifice, by mod- 
ifying a parameter, of Theodorsen's method (referenc-es 17 
and 18) in somewhat arbitrary manner. 

Later it will be shown that the actual pressure dis- 
tribution can be. obtained from the potential theory by 
iteration by first computing the friction layer for the 
potential theoretical pressure distribution and then mod- 
ifying it with the aid of the computed friction layer 
quantities. But, in general, Pinkerton's method in con- 
junction with empirical values, will suffice. 

At small positive and negative angles of attack, how- 
ever, -the pot.ential theoretical pressure distribution 
agrees so clos ely' with .the real to within close proximity 
of the trailing edge that.lt can be used also as a basis. 
At the trailing edge itself the theoretically resultant 
pressure rise toward the . s tagnat i on. point disappears in 
fact (figs, 13 to 19). Theref or e . th^e. pr«s sure distribu- 
tion defining the minimum value .of the profile drag can 
be computed on slender prof ile s . w ithout any measurements 
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Laminar Friction Layer 
a) Pohlhauaen Method 



We employ the known method of Von Zarman and Pohl— 
hausen (reference 2) which has found two more recent and 
critical representations by Prandtl (reference 2) and 
Howarth (reference 21). 

Pohlhausen expresses the velocity profile in the 
, laminar friction layer in the form of a polynomial of the 

fourth rank in ^: ' 

6 

n /n^'' /n\* 

u = a, + - + as ^-y + a3 '^-y ^ a, ^ ^~ . 

where n is distance from the profile and the so— called 
"boundary layer thickness 5 is defined "by 



6' 



J n=6 



= ^o(«) 



(19) 



The fiv«5 coefficients of the polynomial can "be obtained 
from the same number of limiting conditions which are 
§ecured under 
eq^uat i on 



the usual assumption that , Prandtl 's "boundary 



Su Su 
OS on 



i ^1 + V ^ll 



(20) 



where p is air density and. v normal component of the 
velocity in the f rict 1 on . layer , is rigorously satisf.ied 
in points n = 0 and n = 6. They read 



n = 6 : u = u 



0 » 



n = Or u. 



0; ,1^ 



^u 
bn 

0 u 



^u 
d.n2 



> 



— = — u„u 



dn" 



0 0 



(21) 



where ' denotes differentiation with respect t« the arc 
length. The points 0< n< 8 »f the friction layer 
merely satisfy the Karman integral condition of the mo- 
mentum (reference 20, equation 2). 
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ds 

° 0 (22) 



which, simply affords a" statemen.t about average speed val- 
ues. ■ 

If the displacement thickness 6* of the friction 
layer is defined by 

. ". • . h , ^ 

/ / u 

6* = / 1 -dn..- (23) 

/V Uo / 
o 

and the momentum . thickness ^ by 

. = f ~ ^] dn (24) 



O 



the Karma© momentum equation can then be written in the 
form ■ 

-0 - P ~ tuo=-) - »-1f ^ (35)- 

Prom the mathematical point of view it is then ap— 
prop'ri-ate to introduce a parameter A : 

Un 'S ^ 

X (s) = -£ (26) 

V 

which, according to Tollmien (reference 22), i^iay be con- 
sidered as quotient of the local pressure force to local 
friction force, Prom the physical point of view it would 
be better to form it with 5* instead of 5. 

Then the Pohlhausen velocity profile can be written 
in the form; 

■3_ £ x"^— ■■->' 4 - \/.n\^ - 6 - X n\, 

^ ^ '-^ 5 ' y ~ 6"' ' 6 " 2 V 6 /' 2 V S / 6 V 8 > 

■ . . ^27) 

whence follows, if p. d'enotes the air v.iscp.s ity ,. the 
she'aring stress: 
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....\.^.ny 



the displacement thickness 



6 * = 6 



/ 3 



and the momentum thickness 



M. U n(l2 + X) 
6.6 



VlO 120> 



(28) 



(29) 



' A 315 



945 9072/ 



(30) 



S'igure 4 illustrates the velocity profiles for several \ 
values, X = 0 showing, according to (26) since Uq ' = 0, 
the velocity d is tr i.hut i on on the friction sheet, X ■= 7.062 
the velocity profile in proximity of the stagnation point. 
With X = - 12 is afforded, as seen from (28) , the fric- 
tionle.ss. velocity profile charac t er is t io ■ of the "breakaway 
of the laminar frict ion layer. Por X < -12 .the velo-city 
distributions manifest return flow near the wall.. - Alter- 
nately, if X hy great acceleration exceeds the upper 
linri± of Validity -range of the -Pohlhaus en method,. X =.12, 
the result' 'will "be physically irrelevant , "bulged out dis- 
tri"butions (X = 18), 

To single out the physically nonessential boundary 
laye'r thickness 6 occurring in (26) to (30), a common' 
nonlinear differential equation must "be solved, which is 
obtained by inserting (27) into (22). With 



2=1-= 



(31) 



u 



this differential equation reads 



dz 

ds 



f (X) 



+ Z' 



u, 



u. 



g(X) 



(32) 



where 



f(X) = 



7257,6 - 1336.32 X + 37.92 X^ + 0.8 X^ 
213.12 - 5.76 X - X^ 



(33) 
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and g(X) = (34) 

213.12 - 5.76 \ - \2 

The functions f (X), g(X) have teen tabulated "by 
Howarth (reference 21) for the -12 X < +12 validity- 
range of the method. After defining Uq ' and Uq" by- 
graphical or iiumerical differentiation, the isoclinic 

dZ 

field — = const can be plotted and in it the integral 
ds 

curve Z(s), if the stagnation point conditions are ob- 
served, for a given pressure distribution along the pro- 
file and hence a given velocity distribution Ug(s), If 

Z(s) meets the curve Za,(s) = — y at a certain s = 



u 



0 



Sg^, the separation of the laminar friction .layer takes 
place in -this profile point. 

From Z(s) follows, according to (3l), the boundary 
layer thickness 6 (s ) and hence, according to (26), (28) 
(29), (3'^), the characteristic quantities X, t^, 6*, 

and ■» of the laminar boundary layer, A noticeable fact 
is that the diraens ionless . X is a characteristic value. 

for every profile point, independent of Re = ~-- • if 

the pres.sure distribution over the wing contour is inde- 
pendent of Re, as is the case in. the customary test 
and flight conditions. On replacing X for Z in (32) 
a differential equation for X independent of Re re- 
sults. It follov/s particularly that the position rf the 
breakaway (X = -12) of the laminar' friction layer is 
independent of Re. 

It further follows from this irr e l-a t i on that the 
boundary layer thickness and 63 related to two 

different Reynolds numbers Rei * -2ii._iA and Reg = 
are associated through 

U« 1 • 6 ■, ^ 11 - „ ' K 
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By leaving t and . u constant while varying the 
flying speed u^,, 




(37) 



can "be used to reduce the quantities Tq , 6*, of one 

speed for which they were computed with the help of the 
fairly cumhersome isoclinic method to another flying 
speed, 

h ) Note on the Pohlhausen Method ' 

The usefulness of the Pohlhausen approximate method 
is substantially, enhanced "by the fact that it affords, on 
the friction sheet, good agreement with the exact Blasius 
method and on the circular cylinder with Hiemenz ' measure 
ments, More recently, Schubauer (reference 23) measured 
the laminar boundary layer on an elliptic cylinder of 
major and minor axis a^ = 0.299 meters and ag = 0.101- " 
meters, respectively, by means of a hot wire anemometer. 
The direction of the impinging air stream was parallel to 
the major axis, the air speed m^o = 35 meters per second. 

Schubauer obtains from stagnation point — — =0 to 
s .as 

— = 1.832 very good agreement in speed distributions 

a 3 

with Pohlhausen 's figures, but encounters no separation, 

theoretically, for higher values of — — , whereas the ex- 

a 2 

periraent itself indicates separation at — = 1.99. 

Schubauer draws the conclusion that the Pohlhausen method 

is unsuitable by small pressure rise for c orav 'J-t ing the 
separation point when this pressure rise does not last 
long. This inference is, however, inconsequential for 
our laminar friction layer on the wing; in bhe first 
place, the laminar layer will sooner or later turn turbu- 
lent anyhow, before it separates and, in the second place 
the pressure rise las.ts over a much greater distance than 
on the elliptic cylinder. The good agreement of the ve- 
locity profiles in the laminar friction layer with those 
computed by Pohlhausen, on the contrary justifies the use 
of this method on the wing also. 
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The question suggests itself whether the quality of 
the velocity profile could be improved with a polynomial 
of higher tha-n -fourth order .and its coefficient secured 
with the aid of further "boundary conditions. A more ac- ' 
curate computation shows that a polynomial of the sixth 
order gives the shearing stress on the friction sheet to 
within 1/2 percent exactness , that of the fourth order 
only to within 3-i percent. With polynomials higher than 
the sixth order, the coefficients cease to he linear func- 
tions of >.. Because the subsequent boundary conditions 
on the wall are then accompanied by differential quotients 
of the coefficients with respect to s, With every such 
condition containing the derivation of a new coefficient 
in respect to s, the order of the differential equation 
for 5 rises. On top of that, the amount of paper v/.ork 
involved herewith imposes a limit to the use of polynomials 
of higher order which cannot be exceeded even by foregoing 
the compliance of higher boundary conditions on the walls 
and by satisfying only those on the boundary layer limit. 
By these contiguities of higher order on the border of the 
friction layer the velocity profile becomes much more com- 
plete tha,n it actually is. 

Whether the polynomials of higher order reproduce the 
actual velocity profile in every laminar friction layer 
better than those of lower order is itself difficult to 
decide as long as the convergence of this method is not 
proved. 

Attempts might be made to find out whether formulas 
for the velocity profile other than the polynomial in 

~ would give more comprehensive agreement with the exper— 
5 

iment . Kosmodemiansky ' s attempt (reference 24) 



^- = Pf J^sin( 5 J- (38) 



\ 



where ^^^J again a polynomial does not lead any far- 

ther.- 

So for the time being ther-e is no known method for 
computing the laminar boundary layer by selected pressure 
distribution that is definitely superior to. the Pohlhausen 
method. 

One simplification of the latter 's method which may 
be of some advantage in cursory computations might be 
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pointed out : name ly , for a special class of pressure dis- 
t r iTsut i ons the .s 6 lut i on 6(s) , c.an.'be. numerically given 
without plotting df. thfe isoclinic .ft^.ld. If . . 



Uq u^" . 

■ — '■ = C = const ■ (3 9) 

u^'2 



the variables for K in (32) can be separategL, It af- 
fords, according to G-olubew (reference- 25), 



^■'^6. p(X)dX 
Uo + C rT^y 



(40) 



p(X), q(X), ■r(X) "being represented by the rational func- 
tions. • ■ ■ 

p(X) = -213,12 + 5.76 X + X2 ^ 



q(X) = -7257.6 + 1336.32 X - 37.92 X^ _ 0.8 X^ 
r(X) = -213.12 X + 1,92 X2 + 0.2 X^' 



r (41) 



Integration of (3-9) affords the class of the plotted po- 
tential velocity distributions at 

_i_ 

Uq(s) = cs I (1 - C)s + cij for 0^1 (42) 

and 

Uq(s) = Cse®/°^ for C = +1 / (43) 

The integration constants Ci and Cg are defined TDy 
the initial conditions 



^o^^o^ = ^00 and Uo'(so) = u 



1 

00 



For C = 0 and Cj = 0, it affords the stagnation point 
flow; for G = —1., the flow with constant pressure gra- 
dient; for C = +1 and c^ — > oo, the flow on the friction 
layer . . 

A more careful consideration s-hows that it does not 
suffice. to replace a predetermined pressure distribution 

by a polygon 3f straight pieces of — ~. = const (C = — l). 

ds 
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Approximating it by arcs of curves ' on whicK'- -'C • -has '.moften— 
tarily a different value involveB . ju'fi t -as ' mud^ ' paper work 
as the plotting of the- is^ciinic f ield, - S©--- in order t'o 
gain a somewhat clearer insight into the laminar houndary 
layer on the wing the isoclinic field of the Pohlhausen 
method can hardly "be spared at present'.- 

■ The -Trans it ion Point -'• . '■ - ' 

According to Pohlhausen the friction layer on a wing 
can he computed as it would occur if it were laminar up 
to the trailing edge or at least as far as its separation 
point . 

But this assumption is complied with only at small 
Heynolds numbers seldom encountered in practical flight. 
In reality the laminar layer becomes turbulent before 
reaching its separation point. When, as on the friction 
sheet (fig. 3) the transition is assumed to occur abruptly 
in 'a point, the cardinal question of the entire profile 
drag computation arises: namely, where is this transition 
point located. The reason the problem of theoretical 'lo- 
cation of the point is so difficult is due to" the fact 
that it belongs within the ambit of turbulence origin al- 
together. Since, even on the friction sheet' in nonvorti— 
cal flow the transition point defies computation on the 
basis of theoretical considerations, the chances of accom- 
plishing it in an airstream with a certain turbulence and 
accompanied by a pressure gradient over the arc length of 
the body; are even less. Qualitatively it might be ex- 
pected that a high turbulence factor produces an earlier' 
transition becaus.e at least the boundary regions are more 
strongly intermingled and hence the laminar friction layer 
is made turbulent from the outside. It is also suspected 
that, under identical cenditions otherwise, a greater ■ 
pressure rise prumotes an earlier transition'-, since it 
leads to unstable velocity profiles, that is, those with 
ascending return flow risks, and so makes the laminar 
friction layer turbulent from the inside, from the wall. 
Besides, the position of the transitional point is depend- 
ent upon Re = "^^J™* known fro.m the . friction . sheet 
(equation (5))... 

The pressure dist'r i-but'icn over a wing'at constant 
angle of attack in a flow of different air speed is prac- 
tically independent of the ReyholdrS number in. both the 
subcritical and the supejor it ical range. At very small 
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Reynolds numbers the friction layer is laminar- unt il — 
only by greater pressure rise toward the trailing edge — 
it separates in a point, the location of which is inde- 
pendent of Re, If the Reynolds number continues to in- 
crease the friction layer becomes turbulent shortly before 
the separation point is reached at a critica.l Reynolds 

number ^ejj.= -—^ — i which, besides the.wing shape, is 

dependent on the turbulence of the air stream. Transi- 
tion point and laminar separation point coincide. The 
profile drag coefficient drops abruptly because the dead 
air region and hence the pressure drag is reduced by the 
adherence of the turbulent friction layer beyond the lam- 
inar separation point. On further increase in Re the 
transition point travels forward toward pressure minimum. 
That it can advance even before the pressure minimum has 
been only rarely and somewhat uncertainly observed; but 
perhaps merely the Reynolds numbers in the past test and 
flight conditions are still too small to produce such an 
effect . 

For the usual range of Reynolds .numbers therefore, 
it may be said that the transition point lies between 
laminar separation point and pressure minimum. Sharper 
criteria can be expected only on the basis of measure- 
ments at different pressure gradients and turbulence fac- 
tors, which are still lacking for the present in suffi- 
cient quantities. 

G-ruschwitz (reference 26) defined the transition 
point in his plate measurements in the pressure field and 
on the G-Bttingen standard airfoil section no, 387 by plot- 

ting the quotient against- the arc -length, and lo— 

^0 - -' 

cated the transition point where this quotient showed a 

sudden steep rise. On the basis of his experiments he 
found that the Reynolds number 



u^ ^ 

Re^ = ■ (44) 



formed with the laminar momentum thickness * should, for 
the transition point, lie in the region 

360 < Re^, < 6 80 (45) 
Schmidbauer (reference 27) found 
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500. < Ee^y < 790 ' ' " (46 ) 

on c-urved surfaces , and 

■ ■•• ■ • 500 < Ee^- < 550 ' ' (47) 

for the specific case where the transition point is ex^- 
actly located in the laminar separation point. Peters 
(ref eren.ae 28),. ,on, t.h« . other hand., chtained. the- narrow 
int erval 

- ■■ 600 < Ee^ <' 65b' ' " (48) 

f -or- the ' trans it ion' point" on a symmetr ical wing of 16 . 6 
percent thic-knes s . ' Eewr'it ing the transition criterion 
(5) for the flaf plate to Ee^ , a.fferds_ 

.■■R-e^ = .460 (50) 

' y 

since the momentum thickness on. the frict ion sheet is, 
according to Blasius: 




Dryden (reference 29) supposes on the "basis of his ex- 
periments that on the friction sheet the transition 
criterion for ' 0.5 percent turhulence can he raised to 
Ee.,1 = 680' arid for 3 percent turtulence , lowered to 

Ee^ ='210. The criterion (50) is located approximately 
in the center of Dryden 's range. 

Fage 's experiments (reference .30) on Eeynolds num- 
bers in the transition point included circular cylinder, 

airfoil section, sphere, and flat plate. ' He defined the 
transition point as the place nf minimum laminar skin 
friction (fig. 3). This definition is identical with 
that hy Gruschwitis. Page, considering only the hehavior 
of the shearing stress as essential for the character of 
flow and foregoing the determination of the exact form 
of the velocity profile, could- confine his measurements 
t.o the wall adjacent region with linear velocity increase 
rather 'than over -the, total fricticn layer width. .He' 
formed the Eeynolds numher of the friction layer with the 
displacement thickness 6* instead of with the momentum 
thickness 6 and found that 
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* 6 r~~ — 

' -Re|* . u,* =yia , (51) 

was not very dependent upon the tody form and increased 
slightly with decreasing turhulence . Another noteworthy 
fact is that Fage established a definite inflection in 
the pressure d i s t r ihut i ons along the "body contour on the 
transition point, so' that the transition point can he ex- 
perimentally secured without friction layer measurement 
hy a very careful pressure dis tr ihut i on measurement. 

From the wi de ly diver gen t ' range ' of ' Re^ (fig. 5) 
it is apparent that the location of the transition point 
can he given as yet with very little assurance. So for 
the time heing at least the profile drag of a wing at a 

given Reynolds numher Re = — s°_i. should he computed for 

several transition points, which means, compute Re^^Cs) 
with the aid of the laminar momentum thickness defined 
according to Pohlhausen, and select arbitrarily several 
transition points s^ from, say, the Gruschwitz, rjan'ge . 

These transition points are the starting points of the 
turbulent friction layer. 

The Turbulent Friction Layer 

a) G-ruschwitz' Method 

The turbulent friction layer is computed by the ap- 
proximate method proposed by G-ruschwitz (reference 26). 
Like Pohlhausen in his method for laminar friction layer, 
: Gruschwit z starts with the Kafman momentum equation (ref- 
erence 25); but rather than the velocity distribution it- 
self he simply character izes ' it by a form parameter Tj : 



. l--i'^lll 

L ^0 J 



(52) 



One difference existing between the laminar form 
parg.meter X. and the turbulent form parameter ri is 
that X as analytical function' of pressure gradient and 
friction layer measurement is theoretically predictable 
while the calculation of ti dictates other empirical 
laws . 



Since it is necessary to forego an analyti'cal formula 
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for the turbulent velofeity profile-, the'- calculat iQti of 
the four f r ict i on— layer quantities t\ , Tq, 6*, and 
calls for four relations between these four variables 

The f irst represents Karman ' s ' momentum equation 
(25) which , rewritten , reads: . . . 



T. . . d* 1 6 A * ■ d Uo , , 

. = .+ ,f. 1 + - — _ (53) 

p Uq.^ ' .ds . -2 Uq^ ds . 

A second relation .is afforded "by • the .f unct i on 

: -. .H-(n> = r- ■ ■■ ■ " (54) 

•ft 

obtained by Gruschwitz from the evaluation of the record- 
ed turbulent v6loGi-ty prof iles (fig.' 6>. 

The third re la,t i oil was obtaiiied by (JfuscHwitz from 
another; equation deduced from his experimental data. He 
reasoned that the energy change of a'fluid par t icle ,■ mov- 
ing at distance n = 4 from the wall , must certainly be 
dependent on u (^i , s), Uq(s), ^(s), and v. From dimen- 
sional considerations the following form is suggested, 



0 



P 3 

d 



where q = - u. 

. 2 ° 



£ A.| |u^('M + P = - i i- (q n) = Kt). He. ) (55) 
q ds L 2 • J . q ds ^ 

From his measurements Gruschwitz found that the dependence 
of 'Re^ was negligibly small and the' dependence of. ri on 
Ti = 0.8. expressible by a linear law: 

•9 d 

(q ri) = - 0. 00^94 f\ +■ 0.00461 (56) 

q ds 

With the abbreviation ' 

q T, = e ■ (57) 

where-; i indicates the dynamic pressure- difference on the 
friction layer boundary and cn' 6, equation (56) can also 
be written in the form 
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^ 4i = 0.00894 I +-'0- .-00461 a • (5 8) 

ds " ■ ■ ••■ " ' '.' ■• 

The fourth equation is lacking for the time "being; 
hence an estimation of is suhs t ituted. 

With (53), (54). and (58) the quantities n , 6* 

can then TDe appr oximat.ely computed as'frilows: 

Karman's momentum equation (53-) is interpreted as 
differential equation for ^. Posting constant mean 

values for and H = -- and assuming as initial 

P ^0 

momentum thickness of the turbulent friction layer that 
having the laminar friction in the transition point the 
differential equation for * from the transition point 
can he graphically integrated hy a method originating with 
Czuher and described hy Gruschwitz. G-ruschwitz found 

----- =0.002; H = 1.5 (5 9) 

P ^0 

as acceptable approximation. 

The first approximation for ^» is entered in (58) 
for p, which is solved by the same method, by starting 
according to G-rus chw-i t z , in the transition point with n = 
0.1, that is, ,| = 0.1. q. This value of ti in the tran-- 
sition point is arbitrary to a certain extent; above all, 
it corresponds in no manner to the experimental finding 
according to which the turbulence ordinarily starts at high 
n values. G-ruschwitz found, however, that with ti = 0.1 
the theoretical and experimental ri curves agreed best at. 
some distance from the transition point; for the rest, a 
change in initial value does not amount to much because of 
the marked convergence of the curves. According ,t o (57) 
then ri(s) is obtained from r\ (s ) . 

Reading, of f ..the value s of H . for.- .ri- ...f r.om f-igure 6, 
affords according to (54) the displacement thickness 6*. 
This then would leave the shearing stress determination 
in the turbulent boundary layer with thie aid of. a fourth 
relation between quantities 6*, ri , and T^, which 

has not been found to the present time. Gruschwitz 
therefore attempted to define the shear-'ing stress 
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in the following manner: On the ass-umption of the l/ 7 
power law for the velocity distribution in the turbulent 
friction layer 

1/7 

(60) 

V A y 



u„ . 



in -flows without pressure gradient, that is, on the fric- 
tion sheet (uq = 'Uoo ) , the shearing stress can he 
represented, in the form; 

1/4 . 
T / \ 1 /4 

" ° = 0.0225 {---gj = 0.0225 Heg . (&l) 



P -o^ 



Then 



* _ 



8 (S3) 



according to (23), (24), and (60), whence we can also 
write : 

T 



--- = 0.01338 Beg* ^ = 0.01256 Re^ ^^'^ (64) 



p u _ 

G-ruschwitz then proceeded to compute the shearing 
stress Tq(s) in a turbulent friction layer by arbitrary 
pressure gradients, according to (64), where he entered 
Uq(s) and the first approximation of *(s) in t.his equa- 
tion. 

It needs to be proved why and under what restrictive 
conditions this method is justified. ITikuradse (reference 

31) evolved an empirical relation between pressure gradi- 
ent and shearing stress from his measurements in water 
flows with different pressure gradients (wedge flow with 
varying included angle). He found that the speed distri- 
bution could be characterized by a "form parameter" 

1 ds \ V J uo d s S 



Figure 7 .illustrates the relation of the dimens i anles s 
coefficient of the shearing stress 
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■ - J 




1/* 



(66) 



with this dlmeia's i dnl-es's parameter T'l..- It' will "be seen 
that T increases a little at f ir stC w'i'th increas ing 
in. delayed flow (r i > 0), hut then drops rapidly. 
G.ruschwitz' application of the houndary layer law (64) to 
flows with pressure gradients is ohviously therefore 
justified for points on the wall on which Tj departs so 
little from zero that T i (r i ) still closely approaches 
the value ^^^(o) on the f r ic t i on sheet , " In proximity 
of the region of separation (T^ = 0) this inethod there- 
fore f ai Is . 

'In any event, Nikuradee's function T (F ^ ) affords 
a satisfactory substitute for the still lacking, f our t h 
condition between quantities 6*, * , r\ , and To> thus 
enabling the solution of these four factors of the turbu- 
lent friction layer' ofi the wing. 

b) Comparison of the G-ruschwitz' Method with Experiment 

This is 'to' be a brief examination of the extent of 
agreement of the G-ruschwitz approximate solution v/ith 
the measurements on turbulent friction layers. He him- 
self checked his theory first on his own friction layer 
measurements made on the flat plate in the pressure field 
and on a Gottingen airfoil section no. 387; The measure- 
ments the evaluation of which he reported in his article 
were carried out in the range 8.5 x 10 < Ee < 2.6 x 10® 
But he also compared his theory with the Fage and Palkner 
tests (reference 32) on a symmetrical Kar man— Tref f t z pro- 
file of 15 percent thickness, at Re = 1.7 X 10®. The r\ 
curves agreed closely on the whole; but the theoretical 
•fl curve on the Karman-Tr ef f t z profile afforded too high 
values toward the trailing edge, while on the plate the 
theoretical and experimental curves manifested coin- 

cidence in the pressure field. 

St-llper (reference 33) was able -to confirm Gruschwigz 
solution in flight tests at 2.82 X 10®< I^e < 4.88 x 10 , 

and with the same approximate values for — — and H. 

■ -. P ^0^ 

But here also the theoretical i'> curve' toward the trail- 
ing edge seems to lie on the whole above the experimental 
value . 
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Peters (reference 28) exploring the friction layer 
on a symmetrical airfoil of 2.3 meters chord in his check 
of Crruschwitz' theory at larger Ee-ynolds mimhers also 
achieved very good agreement at small positive and nega- 
tive angles of attack, on the "basis of the approximate ' 
valueiB (equation 59): 

— 12~ = 0.0017; H = 1.4 (67) 
P ^0^ 

The lower value of the shearing stress is evidently chosen 
"because of the larger Reynolds nura'ber the amount of which 
is- not given. 

As approximate value for the shearing stress "by a 
given Re the mean shearing stress Tq friction 

sheet at the same Reynolds number Re = — is_ generally 
accepted, 

-JI.^ = 2X (68) 

3 4 
P ^0 • 

c^ being computed from (4) or (6) depending u-pon the chdr- 
acter of flow. 

But Peters' report also discloses the same systematic 
departing of the curves toward the trailing edge by 

increasing angle of- attack, which for the moment cannot 
be explained ; f or the approximate value of T^/p Uq^ ' 
plays precisely in the rear part of the wing where the 
pressure gradient is great, no decisive part; on the other 
hand, the wing itself remains sufficiently flat, so that 
Schmidbauer 's curvature effect (reference 27) can be dis- 
counted. 

A further unusual fact is that G-ruschwitz' H(r)) 
curve was confirmed by Schmidbauer i s well as Peters. 
This is easy to understand since curve H(ti) represents 
a purely geometrical relationship. 



Relations between the Parameter Shape of Laminar 
and Turbulent Velocity Profiles 
This geometric character is seen from the following: 
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In many cases the slieed d=is't^r~iT)ut ieii " in the tiir 

lent ■boundary layer has, while discounting direct wall, 
proximity where a laminar suhlayer with fini.te shtsaring 
stress is formed, appyeximat ely the form of a -power law: 

-f- ^ k = 1/6. 1/7. 1/8 • ■ ■ ■ ■ ■ (69) 



Then function H(ti) must "be given explicit. Because it 
is, according to <23 ) and' (24): . > 



k 

• 8 . " k +""1. 



6 (k + 1) (2k + 1) 



(70) 
(71) 



and 



or 



6 * 

H = — = 1 + 2k (72) 



k = H-=-i (73) 
2 

whence, according to (52), the form parameter: 

(74) 



n = 1 - = 1 



^ 3k - H - 1 



J ■ L H(H + 1) _ 

This function ,^(11); plotted in f igure 6, varies very 
little from the G-ruschwitz experimental curve. 

^ Because of its geometric nature the relation. H(ri) = 
— — is not tied to the flow character in the friction 
* ■ . ■ ■ 

layer and can therefore he applied also tc the laminar 

friction layer. According to the friction layers plotted 
by the Pohlhausen method in figure 6 the parameter ' ti 
varies little from 1 throughout the laminar range and is _ 
therefore not properly suited for representation of the 
laminar speed distrlhut ions . Sruschwitz and Schmidhauer 
"both secured ti values of the order of magnitude of 
Pohlhausen's in the cases where laminar friction layer 
prof iles "were measured. 
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The Point of Separation 



Undei* a strong pressure rise-it;inay happen that the 
turljulent friction Islyer does net adhere as far as the 
trailing edge hut separates hef ore . Gruschwitz adduced 
the separation criterion 



n = 0. 8 



(75) 



while Schmidbauer suspected that" the separation risk is 
postponed to T| ~ 0 . 9 on the assumption that the pres- 
sure rise at high f| is flat. Peters' findings comply- 
very closely with the G-ruschwitz separation criterion 
(75), if the r\ values secured from the measured speed 
r! is tr ihut i ons are used as a "basis. The theoretical solu- 
tion of ri , however, affords, infer ior values in prox- 
imity of the separation pointj hence no separation is 
at all urediotah le therefrom. 



Gruschwit? likewise associated the location of the 
separation x)oint with that of the transition point. The 
earlier- the transition takes place the greater becomes 
the turhulent houndary layer thickness in a fixed profile 
point hehind the transition point. And this is accom— 
panierl hy a greater risk of separation; hence at constant 

Re = — — — the separation point t) = 0.8 movfes forward 

when the transition point ig shifted forward. By fixed 
transition point and ascending Reynolds numher, on the 
other hand, the separation region shifts toward the trail- 
ing edge because the momentum thickness 'becomes less. 
These two effects are superimposed on the wing in such a 
way that the effect of the Reynolds number on' the separa'- 
tion point follows as the difference of two effects, the 
Torefix of which does not appear to 'be a'bsolute. At small 
and medium Reynolds numbers the effect of transition point 
travel usually -or eponderat e s , that is, with ascending 
Reynolds number the transition •ooint shifts forward, the 
drag increases ; But , when at larger Reynolds num'bers the 
transitional region is already far forward so that its 
location can -be looked .upon as appr oximat e.ly c onstant , the 
first effect dis-appears- and the separation, point. shifts 
backward with ascending Reynolds number., the .drag decreases. 

TTikuradse . (see fig. 7) also secured a separation cri- 
terion from his experiments, at least by extrapolation: 



r a = 0.08 



(76) 
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which is in go od agr e ement> w i t h the (Jruschwi.tz criterion, 
as will "be , shown,, later . , . 

Since the' dynamical effect' of the dead air region is 
to he disregarded, we possess- herewith the theoretical 
raean.s to- compute the friction layer f-rdm stagnati-^on point 
to separation point or on thin profiles as far as the 
trailing edge. 

IV. METHOD OF PREDICTING THE PROFILE DRAG PROM 
THE CHAEAOTBEISTIC.S OF THE FRICTION LAYER 
1. Friction Drag 



The friction drag is determined according to (l) hy 
the variation of the shearing stress Tq, which itself 
is a characteristic quantity of the friction lay-er . 

Karman's momentum equation (25) represents a general 
relation "between Tq, 6*,. and If the two character- 

istic lengths 6* and ,i could he computed without re- 
course to the Karman equation, this would afford- the pos- 
sibility of computing the shearing stress and so, the 
friction drag. Integration of the shearing stress result- 
ing from the Karman momentum equation along the profile 
surface gives the friction drag per unit length of span: 

^r = y T 0 cos ^ d 8 =. p[uo^ ^ cos ^] j^jj ~ S ^ ^ ^ ^ 

I+II I+II 

(7?) 

The suhscripts under the integral signs, again, denoting in- 
tegration from stagnation point with respect'to su<j.tion 
side and pressure side, and the suhscripts on the square 
brackets the corresponding values of the "bracketed trail—, 
ing edge term. 

In the laminar friction layer the shearing stress can 
he approximated by the Karman— Pdhlhausen: method. But in 
the turbulent layer where only equations (53 ) , (54) , and 
(58) are available for defining S*, ,^,11. and Tq , an 
estimation of the shearing stress jtj has to" serve as 

fourth equation. For the solution of and 5* the 

unc-ertainty of T plays no great part, since the'se quan— 
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titles geiierariy are not very' dependent on . Tq,' 'bTit- in 
the determination of the friction drag- the' uncertainty' 
of the Tp appraisal is significant according to equa- 
tion (64). A r« liable solutio.h of T.^ . and hence of the 
profile drag could be. secured if • t-^ could be t-ied in 

with the form paramfeter ri , as on th'e laminar friction 
layer . . • " ■ 

2. Pressure Drag 

The pressure drag of a wing has been determined up 
to now by plotting — according to (2) — ' the pressure p 
against the coordinate z of the profile pointsnormal to 
the direction of flight and defining the area of this 
pressure distribution. However, as has been pointed out 
by Betz (reference 6), the theoretical' uncertainty of this 
method is that this area consists of comparatively large 
positive and negative parts, as a result of which the 
pressure drag is obtained as difference of two doubtful 
values of the same order of magnitude. This defect has 
been > subs eq.uent ly substantiated from various sources (ref- 
erences 34 and 35). 

A method is therefore developed which is amenable to 
the theoretical solution from the friction layer quantities 
and obviates this uncertainty. Strictly speaking, v;e con- 
fine ourselves to the case vhere no separ-atioa occurs in 
the friction layer; hence the, pressure drag is exclusively 
due to the fact that the adhering friction layer pushes the 
potential flow away from the surface to the amount of the 
displacement thickness 5*. But if s eparat i on ■ occur s , the 
effects of. these separated air balls on the wing must be 
taken into consideration, according to Eetz (36), As long 
as the separation is insignificant, these reactions are 
minor and the pressure drag can then still be estimated, 
at least approximately. 

The line of reasoning of this method ig as follows . 
(fig.' 8): ■ 

Visualize a flow past the s,urfac.e which is friction- 
less, but with the same displacement of streamlines in 
relation to potential' flow as occurs on the natural fric- 
tion—encumbered flow. This displaceme.nt is- effected by 
superposing a source on the surface. With this superposi- 
tion of source flow on the original potential flow a new 
potential flow with the same pressure distribution as the 
actual 'flow is obtained. 
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The control area on which the force 'eq.uiliTDriiim is 
to he analyzed, Is then defined as follows; Place a 
skin & at equal distance 6* around the wing T, 5* 
itself heing defined hy (23). This skin represents, as 
it were, the surface of the displacement water through 
which no fluid of the source flow passes any longer. On 
the trailing edge this skin is closed off hy a plane H 
at right angle to the direction of flight. The pressure 
and momentum on the closed control area formed hy T , G- , 
and H are now analyzed. With D denoting the resultant 
in flight direction stemming from the pressure forces and 
that produced hy the momentum with I, mu"st follow: 

+ Ijv. .+ D(, + I(, + + Ig = 0 (78) 

Area G- having "been chosen so that no momentum passes 
through, it is therefore 

I(j = 0 (79) 
whence the pressure drag of the wing follows at: 

-Dj = + + Ijj + Ij, (80) 



The solution of the momentum entering through the 
v;ing surface and plane H at the tra iling ' edge must he 
preceded hy the determination of the source superposition 
on the wing surface. At a selected point s = s« of the 
wing profile the fluid volume UqS 
length of span per unit time. At 

d(uo6*) 



is displaced per unit 
ds this vol— 



s = s 



ume amounts to 

hy the friction 

at s=SQ+ds 
terpreted as if 



u^5* 



ds 



ds . Hence the displacement 

layer per unit time is 

This fact can he in— 
span direction a source 



greater hy ^^i-l 
ds 



than at point 
on line s = s 



0 • 

in 



had heen superposed with a yield per unit length of span of 



dE = 1^1:1 ds 
ds 



(81) 



Hence the momentum per unit time enter ing at s 
region hounded hy T. , G- , and H, amount.s to 



Sg, th^ 



dig,' 



d(u-S*) 

^ — p Uq d E = - ^ Uq — : ds 

ds 



(82J 
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where entering, impulses are negative, transported impulses 
are positive. The total impulse passing through the wing 
surface is . 



It. • -P / — 



■) 



is 



(83) 



i+ri 



and its component in flight direction 



Irj, = - p 



till 



O P 

ds 



r d(u^6*) 
^ J ° ds 



dx (84) 



1+ I-I 



Through the plane at the trailing edge .the momentum 



is transported, vhere Uq , Uq 

Hi Hi I 

locities at the trailing edge on suction and pressure side 
and Ej, Ejj. the corresponding source yields: 



(85) 

are the potential ve— 



p d(u^6*) 

/ ds = u 

J ds 



E 



II 



II 



d(uQ5*) 
ds 



0 Ht** Hi 



ds = u^ p 8*rr 



(86) 



Its component in flight direction is 

= P C^o H2 5h* cos Ph 3. 



I+II 



(87) 



with * ^Hjj denoting the angle hetween profile tan- 

gent and flight direction at the trailing edge. 

In order that the -pressure drag may he computed, 

the pressure force on area G- and plane H must he de- 
fined. The first follows from Bernoulli's equation at 



I+II 



(88) 
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. The second from the argument: The continuation of 
area 0- . placed at distance 5* around the wing, 'beyond 
the trailing edge as area S so that it's width is every- 
where equal to the displacement thickness 6* of the 
dead— air region, does- not alter the potential flow out- 
side this- surface nor, hence, the pres sure . di s trihut ion . 
However, G and S represent a so— called conoid on 
which the resultant pressure force of a potential 



is 
■or 



known to' Ve' "zero ; that is, it affords'. 



+ D, 



flow 
= 0 



(89) 



With 35* (see footjaote) denoting the width 
hody at any point and p the pressure, the 
force in the wake flow becomes 



of the wake 
pr es sure 



6c» 



= 2 



p d 6 



(90) 



where 2 



is the displacement thickness for Uq = u 



Therefore the di s-olacement thickness in the wake 
with respect to the outside pressure must he established. 
It ca,n be assumed approximately that over- the short- 
distance in which the pressure abates to its final value: 
no ap-oreciable energy changes take place in the wake, 
so that each particle is merely subject to the accelera- 
tion due to the pressure differences. Then to each 
stream filament Bernoulli's .equation 



^H^^ = I ^^0^ - ^0 = Ph - P 



(91) 



u(n) 



can be applied and hence the wake configuration 
at every other place in relation to the prevailing pres- 
sure computed from the velocity profile at the trailing 
edge ^gdijj). 



The velocity profile in the friction layer at the 
trailing edge might be approximated by the power law 

= ^0 h(^0''' (92) 

It is 25* = 8 J *+ 5 J J*, where 6j*, 6jj* indicate the 

proportions of the. wake displacement thickness, emanating 

from the friction layer flows at suction and pressure 
sit^e; on the trailing edge itself 26 g* = 6 g' * + 5 g *. 
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To simplify- mat t er's ; ■ i t' is a-ssu-me"d' that" tlie flow is 
almost" symmetrical-, "that , the 'saii^ on pr fes sure-' s ide 
arid suction si'de. ' The exponent ' k arid the hcundafy layer 
thickness 6jj f-ollo-w from the displ^ceme'nt ' thickness 

aiid the .momentum thickness ^ on tie trailing edge . ac- 
cording t c . (70) arid. (73 ) 

Equations (91) and (92) aifford first the wake veloc- 
ity on a stream filament which on the trailing edge is 
hy ng distant from the wing surface. 



By given Ug it affords in the center of the wake 
hody (n = njj = O) : 



u — 

^^0 H 



^ ^- - 1 (94) 



The general coordinati on of n and ^jj » which first 
defines the vel oc i ty pr of ile u(n) and then the displace- 
ment thickness 6* of the wake, is found in the continuity 
equation: 

- -u d n = u„' d' ri„ (95) 
xl n 

when starting from n = 0, we proceed- step >y step from 
one streamline to the next. 

For k = 1 and k = l/ 2 , " u(n) can even he obtained 
hy sq.uaring. With the ahhr e viat i en 

---S- = Y • (96) 

U TT 

0 H 

equation (95) gives for k = 1: 



u 



^o H ■ 



^ +/y^ - 1 (97) 
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and for k = l/s 
^0 H y ^o' H 



n 

8h 



V 1 + i 



- (Y^ - 1) In 



1 

r.*-.^ 1 




~ — ^ Y^. + 1 + 
o H 



(99) 



The curve of half the wake "body width is shown 

plotted against ■ Y' 'for. fc = 1, l/2, l/4.,' and. if 8 in 
figure 9. ' . • • " ' . . • , . ■ 



Furthermore, with the pressure in the section of 
width 2 6 * at 



(100) 



the pressure drag !)<; exerted en the wake follows from 
(90) at: 



/ 

6 * 
H 



p.d 5*,= p u^ 6jj* Qi{ k. 



(IQI) 



Factor Q, is illustrated in figure 10. Exponent k cor- 
responds to G-rus chwrtz ' parameter r\ "because of (74), 
according to the tahle: 



k . 


. 1 ■ 


1/2. 


1/4 


1/8 


n 


0. 973 


0. 833 


0.633 


0.454 



Hence, with the integral in Irp partially integrated and 

angle ,$ disregarded', ecLuations (80), (84), (87), (88), 
(89), and (lOl) give thfe ' pres sure drag at " ; 



J \ ds 
I+II \ 



L • v^o H y • 



- Q) 



(102) 
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or, with the a'bbreviation, 

= M; ("k. =f _!i«_y (1 ^ _ 1 . (103) 

finally: 

I+II 



Herewith the calculaticn of the pressure drag has teen re- 
duced to that of the displacement thickness 6*, that is, 
to a friction layer quantity* According to figure 11, fac- 
tor is little dependent on the shape of the velocity 
profile at the trailing . edge , "but more on' the ratio 

^0 H 

Addition of (77) and (104), with angle 3 disregard- 
ed, approximates the profile drag Wp of the wing to 

% ~ p u^ jj (2 * g - \t; 5g*) (105) 

Equation (105) permits particularly the reduction of 
the profile drag measurement in practice to a friction 
layer measurement on the trailing edge and therefore rep- 
resents a welcome supplement to the Betz momentum method 
(reference 37) . 

But for purely theoretical solution of the total 
profile drag, equation (105) is still inapplicable for the 
present, "because its principal constituent is the friction 
drag, and this will continue until the momentum thickness 
4 in the turbulent part of the friction layer can be de- 
termined without approximate assumptions on the shearing 
stress Tjj. 

So, for the time being, the procedure of defining the 
profile drag will have to consist of computing as 
well as 5* along the. profile tlyontour with the -methods 
described in chapter III and then of the friction drag 
with (l) and the pressure drag with (104). 

Before proceeding to the practical application of 
the proposed method, v/e v/ish to point out how the method 
may be employed to solve the real pressure distr iliut ion 
from a given potential pressure distribution, by itera- 
tion. 
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Start ing f r'om tlie. . ;^otgat ia;3t pa; $-\irer dis tr ibut ion , 
we compute with it t he dis.J)lacement.-,t-h.ickness • • "6 * ; the 
curve ' having the distance 5*- "-from the initial profile 
represents a new profile c ont our " whos e potential pressure 
distrihution is obtained hy transformatian- on the cir.ole 
and presents a first .approach to the-ac:tuai pressure dis- 
tribution of the ijttitial profile. ,rhe.n the displacemeht 
thickness .8* -for this pressure distrib.ut ion is computed, 
affording a new profile contfiur, etc. 

Now the developed method is being used to define the 
profile drag coefficient on the seven symmetrical Karman— 
Trefftz profiles studied by Page, Falfcneu , and Walker 
(reference 35) in symmetrical flow. 

V. APPLICATION TO THE SOLUTION Of PEOFILE DRA& OP SEVEN 
SYMMETRICAL KAEMAN-TREPPTZ PROPILBS IN SYMIIETRICAL PLOW 



These profiles were chosen because their profile .drag 
and pressube distribution was recorded at 6.58, 10.76, 
15. ."^2, 19.9, and 24.25 meters per second air speed; besides 
friction layer data for one were available (reference 32). 

Thickness and trailing edge were so defined as to 
locate the maximum thickness at l/S wing chord. Then 

0a."02 

K-= (106) 

... 2a.: 

according to Grlauert, but .;in ,.Betz' method of identifica— . 
tion (reference .16 ) , where k defines the trailing edge 

cp = K n (107) 

and , . . • ■ 

■ ' .0 0 ■ - .. ■ . . . > 

■•€ = .--i— S =.2 K . . , .(108) 

a ■ 

the thickness. The actual chord of the models in the 

0 ,2 • 
tests was smaller by — — — percent than the therretical 

€ 

chord; the trai ling edge was r ounded ■ of f . The form param- 
eters are given in table 1; the profiles themselves (fig. 
12) are replicas of those -in reference (35).; 
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■ TABLE...1.- PE0FILE PARAMETERS 



Profile 


e 


K. ,. 


cp 

(deg) 


* the or 
(m) 


^ac ttia 1 
.(m)' - 


(m)- 


d 


1 
2 
3 
4 
5 

6 ■ 
7 


0 . 0333 
. 0667 
. 1 

. 1429 
.2 
.25 
.333- 


0 . 0167 
. 03 33 
. 05 
. 0714 
. 1 

. 125 
. 1667 


3 , 

6 

9 

12 . 9 

18 

22.. 5 
30 


1.462 

1.51 

1 . 03 
.7475 
,564 
.47 
.373 


1 .375 
1.465 
1.009 
. 736 
.559 
.466 
.371 


0. 0756 
. 1524 
. 152 
. 1524 
. 1524 

.1524 
. 1492 


0.05.51 
. . 1040 
. 1506 
.2069 
.2726 
.3270 
.4025 



is 

The kinematic viscosity of air/ computed from 
English data on Re, Uoo, and t, at u = 1.48 X 
square ' meters 'pe r second. 



the 
10-^ 



From' the pressure distributions plotted in figures 
13 to 19, it is seen that, apart from profile 2, thedis— 
tributions vary so little and with ascending Reynolds num 
"ber so uhsys temat ically from one another that the theoret 
ical distribution, even on thick strut sections, may he 
lookedupon as satisfactory up to 80 percent of the chord 



With the aid of the theoretical pressure distri- 
bution the laminar friction layer was computed for the 
lowest wind speed of vl^ = 6.58 meters per second ac- 
cording to the Pohlhausen method. As direct graphical 
result of the isocline method, figure 20 shows the val- 
Z u_„ /R \ ^ 

Ties of = ( - i Re (see (31)) plotted against the 

t \^ty 

arc length up to the value where the theory stipulates 
separation in laminar friction layer. In figure 21, the 

location of the separation point cf the laminar 

layer is plotted against profile thickness d, along 
with the correspcnding values for symmetrical Joukowski 
profiles according to Zosmodemiansky (reference 24) and 
the separation point for the circular cylinder (d = l) 
by Hiemenz / ag \ 

and the circular cylinder ( d = — = 0.338 by Schubauer 

• V ^1 / 

Prom Z,. the Pohlhausen boundary layer thickness 
6 can be obtained with (31), so that the charac-ter is- 
tic quantities of the laminar fr ict ion. layer X, 6*, 
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and Tq (figs. 22' to 24)' can be c biput ed f r om' (26 ) ., (28), 
(29 ) , and <30) ; The corresponding values of 5*, *, and 
T^"for Tiw =^ '15 .SS fiieters per'second and Ucs* 24145 

meters .per, second y.iji'd' speed , as af f or ded ' f r bm (3?) "by 
conversion, from the values for ' Uco = 6,.5p, meters' per .sec— 
ond'are shown individually ' in f igures" . 25 1 6. 36 .. 

For .the- select.ion of the trans it ic.n. points the 

Reynolds num'ber Re^ = -~- — was plotted against the. arc 

length with' speed' .''Ugc^. .as parameter, and tRen. several 
points chosen arbitrarily .from G-ruschwit z ! • 4''iO't ed Ee^ = 
range . '.■•.,•• 

Figures 25- to 36 illustrate the momentum thickness 
the form parameter n . and the displacement thickness 
6* in the tur'bTiilent friction layer conformal to G-rusch— 
witz' solution for profiles 1, 3, 5, and 7 and -air speeds 
uoo = 6i58, 15.32, 24.45 meters per second. at different 

transition point locations. Before the transition points 
are the values of the laminar friction layer obtained by 
Pohlhausen theory. 

That * and 6* in the transiticn point are com- 
puted identically is due to the fact that the basic val- 
ue Ti = 0.1 in the transition point defines a value of 

H = — — ~ 1. It is known from experiment that the dis- 
placement thickness on the transition point frequently 
grows very rapidly almost discontinuous. Since the amount 
of this discontinuity cannot be c omput ed, quit e apart from 
the fact that the T) values immediately behind the tran- 
sition point are themselves uncertain, the curve of 5* 
over a certain distance behind the transition point is 
not accurately determinable. For the calculation of the 
pressure drag t,he laminar displacement thickness was 
therefore continued with a somewhat arbitrary curve up 
to the more safely ddfined 6* values in the rear por- 
tion of the profiles i These connecting curves are not 
shown . 

As 'the o.harapt.er.i-.st ic:. -quan't it ies- in the .turbulent 

• ■ - ■ Uco t ■ ■ • 
friction layer cannot • be. reduced from one >B.e = —- — . 

to another, t.hey must be graphically, and numerically de- 
fined for each Reynolds number. 
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The momentum thickness in the turhulent layer 

is seen to grow so much more at constant w.ind. speed as 
the separation occurs earlier. It is also apparent that 
the form parameter- :T) at. some dis tance from ' the transi-* 
tion point is not suhs tant ial ly affected "by a displace- 
ment of the transition point -nor •'by the Eeynolds numher. 
It may therefore 'he -asserted that, in the range of Ee 
considered here at leas.t , the region of separation of the 
turhulerit friction layer on the wing is itself not very 
profoundly affected hy the Reynolds numher' (fig. 37), as 
suspected, in fact, hy Lyon (reference 38); it is li.ke — 
wise apparent- from the' experimental pressure' distrihu— 
tions whose departure from the theoretical pres'sure dis— 
trihution afforded a first rough estimation of th-e loca- 
tion of the separation point. 

Jigure 37 ftirther discloses the f orward " s hif t of- 
the region of separation ' with growing profile thick- 

ness , 

An unusual fact is that the separation criteria of 
Oruschwit 2 - (75 ) and Nikuradse (76) yield the separation 
point at the same place (fig. 38). 

In figures 39 to 50 the values of the shearing stress 

■^0- X 

-— are shown for three wind speeds plotted against — 

" t 

P Uoo 

for laminar and turbulent friction layer. The values of 
the turbulent friction layer were throughout computed 
according to (64). To illustrate the manner in which the 
sh-ear ing 'st'res s diverges from that of the frict.i.on sheet 
•it was als'o jplotted for the same Seynolds number Ee = 
"Uqo ■ t " ■ 

. This estimation yields, on the whole, too low 

V ■ 

drag values . 

The shearing stress distribution on profile 3 at 
u^= 34.45 meters per second wind speed (fig. 48) reveals 
a good agreement of the theoretical values of t with 
those obtained from the friction layer measurements (ref- 
erence 32 ) . 

Prom the shearing stress curves, which at greater 
profile thickness must-be modified so as to make this 
stress ' disappear in the separation point of the turbulent 
friction layer, the friction drag is obtained by itera- 
tion according to (l). Because the drag portions of 
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s'u'cti on" and pressure 3 ide are equal as' a result of flow 
symmetry j the . f r. ic t i on drag coefficient reads:- 



• CO 



In figures 51 to 5*? this d bef f icient - o^v " is. shown 



r 



plotted:.-against .the coordinate — — of the transition 

■ t ■ 

point for all seven profiles at varying wind speeds, 
along with the experimental fription drag factors of 
Fage , Falkner, and Walker after suTstracting the pressure 
drag from the weighed total drag at Cg^ = 0. Since the 
transition point of these tests is not known, except for 
profile 3, its location was deduced conversely hy plot- 
ting the experimental c^^^ value on the theoretical 

curve for ''wr^^~y • wind speed increases the 

transition point on all profiles moves forward, but with- 
out overstepping the pressure minimum. On profile 3 and 
at Uoo = 24.45 meters per second the experimentally de- 
fined transition point is located at Xu./* = 0.24 com— 

pared with -- = 0.27 for the experimental c„ value 

in figure 53. Figures 51 to 57 further contain the co- 
efficients of the equivalent friction sheets according 
to Betz' and Jones' approximation (sec. II) where ^xJ^- 
is to he computed hy (5). Comparison discloses that 
Beta' equation (7) already gives a very good replica of 
the friction drag for medium wing thickness.. 

The pressure drag coefficient fallows from (104) at 
o 

A careful study discloses that the pressure drag 
coefficient also drops "by constant Reynolds humher if 
the transition point is. left to shift "backward, although 
this relation is too insignificant to, warrant graphical 
representation. Figure 58 shows the theoretical and ex- 
perimental pressure drag plotted against profile thick- 
ness at three wind speeds; the agreement is satisfactory 
as long as the friction layer does not "become separated. 
The consistently slightly lower theoretical pressure drag 
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values are likely to "b'e due t-Q the fact ■ that- the G-rusch- 
witz solution yields unusually, high moment and- displace- 
ment thickness values toward the trailing edge, as has 
been pointed out. Figure 59-shows the change in pressure 
and friction drag "by constant Reynolds number, as copied 
from reference (35), 

In conclusion, it can he stated that with the avail- 
able data on the flow in the friction layer, the profile 
drag, can- be approximately computed, provided the location 
of the transition point is accurately known, 

VI, GEHBEAL CONCLUSIONS 



As shown in section V, the location of the transi- 
tion, point is of decisive- importance for the theoretical 
solution of the profile drag. Since the available data 
for a positive identification of this point are inade— ■ 
q.uate , the available experimental findings in conjunction 
with the results of the present calculation are employed 
in an attempt to find this point even. if only very rough- 
ly.. 

The suspicion that the location of the transitirin 
point is related with the pressure distribution over the 
surface of the body suggests the study of relationship 

of Reynolds number Re.^ = """^J^ ^'■^^ transition point 

with the dimens i oniess quantity -— — Ji2 in the transi— 

"0 d s . 

tion point, SruschWitz (reference ' 26 ) labored on. a so- 
lution of this problem but was unable to establish a re- 
lation between the two quantities, J"age .(reference 30) 
also merely found that the Reynolds number Re u, = 

0 

i * ■ 

-— in the transition point is not very much depend— 

ent on the pressure distribution. 

Therefore, assuming that the compariBcn of the ex- 
perimental friction drag factors on the tliin Karman— ., 
Trefft.z profiles of Page, Palkner , and Walker (reference 
35) with our theoretically computed values , gives the lo- 
cation of the transition point confdrraal' to figures 51 
to 57 sufficiently- accurate affords' additional data for 
exploring the relation of Re.^ in the trans it i cn point 




and' -- ^ '" in the transition point. The satisfactory 
Uq 4 s 
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'"accuracy '"^6 f our calciilat'lon 'may -¥e • consider proved "by 

the fact that on profile 's the experimentally defined 

transition point (reference 32-)' I'ie's at' = 0.24' t 

against . our = 0.27 t i' 

Figure 60 illustrates' how on 'prof iles '1 to 6 at 
varying wind speed- the dimens ionless quantity 

4 d Uq 

— — — changes with Re. over the profile contour, 

u-o d s ■ . . 

the intersection points of the curves with the Re^ 
axis def.ihing the pressure minima on the profile. iEach 
curve shows the point in which, according to figures 51 
to 57, the transition occurs. Also shr>wn are -the corre- 
sponding values in the transition points o'btained "by 
Gruschwitz and his test series IV, V and airfoil section 
no. S87 at a = 12°, those "by Stiiper (reference 33) for 
his test series I to 17, and those "by Page and Palkner 
(reference 32) on profile no. 3. 

It is evident that this mul t ipl i c i t'y of points es— 
tahlishes no definite relations-hip "between Re ^ and- 

"i "iio 

— in the transition point. But it is startling 

Uq d s 

that the pressure gradient is apparently quite unimpor- 
tant as far as transition is concerned and that the value 
of is largely decisive. All the Re values of 

the transition points are within Dryden's range of 210 < 
Re^ < 680 (fig. 5). The reason the "street" of the Re^ 

values is so wide may have its "basis in the different 
degree of tur"bulerici3 of the air stream and the different 
surface roughness. Small Re^ values in the transition 
point signify high turhuience or c ons id.era'ble roughness, 
large Re^ values, low tur'bulence and pmooth surface. 

If our pr esumpt i on that, the transition point is in- 
dependent of the pressure gradient and merely .affected "by 
tur'bulence and roughness is' correct, the.n the transition 
point must also "be a"ble to shift "before the pressure mini- 
mum h.y constant. tur"bulence and given roughness, 'when the 

Reynolds numher Re = --— - is large' enough. For, from 

(36) follows that 
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in a fixed profile point. Jrom the test data availahle 
this fact has "been -.scarcely recognized whence the general 
opinion that, the transition point did not shift hey'ond 
the pressure minimum. The present analysis at any rate 
makes it plausible that the pressure minimum does not 
form such a boundary for the transition point and_ that 
the reason such a location cculd not be satisfactorily 
ascertained is solely due to the too low Eeynolds numbers 

Re = — — — — of the experimental material. 

u . 

A test at Re = lO''' made in the meantime appears 
to confirm this assumption. Gruschwil; z ' measurements 
ranged between 8.5 x 10^ ■ < Re <2.'6 x 10^, Stiiper's. 
between 2.8 X 10^ < Ee < 4.9 X 10^, and Page, Palkner , 
and Walker's between 1.7 x 10^ < Ee < 2.4 x 10^, 

Figure -61 shows Ee.^^ pl&tted against the arc length 
of profile 4 for five different Reynolds numbers. The 
lowest three curves represent the Fage , Falkner , and 
Walker test series at ' u =. 6.58 meters per second, 15.32 
meters per second, and 24,45 meters per second. 

The transition points obtained with the criterion 
(50) He^5 = 460 are located "behind the pressure minimum. 
But by a tenfold or a hundredfold increase of the lowest 
Ee the criterion Ee,, = 460 would bring the transition 
points before the pressure minimum. 

In conclusion, our presumption that the transition 
point is simply characterized by a certain Re,n value 
which itself is dependent upon turbulence and roughness 
is to be used to compute theoretically the behavior of 
the profile drag at much larger Eeynolds niimbers than 
ever reached in the test. Base the transition criterion 
by moderate turbulence on the condition of, say, Ee = 
460; then the transition point, and hence the turbulent 
friction layer, can be computed with the aid of the 
laminar momentum thickness. The calculation, involving 
the profiles 1 to 4, consists in ascertaining the two 
portions of the profile drag at Eeynolds number exactly 
l/lO or 2/10 percent larger than the lowest experimental 
Reynolds numbers ^e^^^ used by Page, Palkner , and 

Walker. Since at this divergent Eeynolds number G-rusch— 
witz' approximate assumption (59) for Tq can no longer 

be used to compute the turbulent friction layer, the new 
average values 
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Re 



E.e_ 
min 



10 



100 



P ^o' 



0.002 0.00175 



0 . 00125 



are comp-uted according to eq.uations (6) and (68). 

It appears that at large Reynolds numbers the pres- 
sure drag coefficient drops "but very little with increas- 
ing Re, while the friction drag coefficient continues 
to decrease uniformly. Because of it the percent portion 
of the pressure drag would rise c ons idera"b ly with ascend- 
ing Re. This shift of the portion of friction and' pres — • 
sure drag with increasing Re would also favor the most 
appropriate shape of the wing to the extent that with as- 
cending Re a decrease in pressure drag becomes increas- 
ingly more important, that is, more slender profiles 
would have to be used with increasing Reynolds number. 

In order to be able to remove the uncertainties .still 
afflicting the present calculations the principles of the 
theory require a still better explanation in some points; 
the most urgent problems are; 

1. A more definite identification of the location 

of the transition point 

2. A more reliable arithmetic method for solving the 

shearing stress in the turbulent friction 
layer 



3. Clarification of the discrepancy between experi- 
mentally and theoretically defined momentum 
thickness in the turbulent friction layer near 
the trailing edge. 

Not until these problems have been solved satisfac- 
torily will it be possible to compute the profile drag 
of a wing in purely theoretical manner with the necessary 
accuracy . 

The drag proportions for the four profiles then are: 
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Jigures 1,2,12,13,14,15,16,17. 
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Figure 13. -Pressure distributions. 
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Figures 3,4,5,6. 
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Figure 4. -Velocity profiles in 

laminar friction 
layer (according to Pohlhansen) 



Iv 



H- 



Pohlhauson 



/ 



^'G-ruschwitz 



-* TV* 



TTTo" 



T,0^ 
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on turbulent velocity profiles. 
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Figures 7,8,9,10,11. 
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rigure 18-19. -Pressure distributions. 
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Figure 20.-GrapMoal result of 
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method. 
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Tigore 22.-IiaminAr form parameter X. Figure 24. -Shearing stress in 

laminar friction layer. 




Figore 25-28.-6*,^ and n in laminar and turbulent friction layer for 
profile 1,3,5 and 7 at uqo = 6»58 m/a. 




Figure 29-32.-8*,^ and 11 in launinar and turtulent friction layer for 
profile 1,3,5 and 7 at = 15,32 m/s. 




Jigare 33-36.-6*,^ and r\ in luminar and txirbulent friction layer for 
profile 1,3,5 and 7 at = 24.45 m/s. 
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Figur* 39-42. -Shearing stress tq in laminar and turlwlent friction layer 
for profile 1,3,5 and 7 at 6.58 m/ %. 




Vigor* 43-46. -Shearing stress Tq in laminur and turbiilent friction layer 
for profile 1,3,5, and 7 at Ugj* 15,32 m/s. . 
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Figures 51 .52, 53 , 54, 55, 56 . 57. 
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Figs. 59,61 
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Figure 59. -Profile drag coefficient plotted against 
thickness at two Re. 
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